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' Abstract 

c ; 

O . Details of the dynamic behaviour of different charged clusters in 

liquids are discussed. Their associated mass is considered which pos- 
sesses a number of interesting features in a normal viscous liquid. 

> 

I In the two previous papers [1, 2] the authors noted that the systemat- 

ically observed [3-6] temperature dependence of the positive ions effective 
^ ' mass in superfiuid helium has mainly normal (non-superfluid) origin. For- 

Q ■ mally, the point is that under non-stationary conditions the Stokes drag force 

^ . acting on a sphere moving in a normal fluid actually has both real and 

Q I imaginary components. In the extreme case of low Reynolds numbers one 

has [7, 8] 

>< : F{u;)= GTir^R (^1 + ^) v (u) + 



where p is the liquid density, r] is its viscosity, v is the sphere velocity, 6 (u) 
is the so-called viscous penetration depth, and R is the sphere radius (or 
effective radius of the ion or a different particle or cluster). 
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It is natural to identify the coefficient at iu!v{u!) witfi the efficient associ- 
ated mass of the cluster: 



m''''{uj,R)=mid{p,R) 



' 9 5(cu) 
1+ ^ ' 



2 R 



m,d{p, R) = 27rpi?V3 (2) 

The associated mass m"** proves to be frequency-dependent (the dependence 
being rather strong at low frequencies) and this circumstance should be ac- 
tually taken into account when considering the ion (cluster) dynamics em- 
ploying the Navier-Stokes equation. 

Experiments [3-6] were carried out at finite frequencies, and the data of 
Ref. [3] were obtained in the range of ujt < 1, where r is the the particle 
velocity relaxation time. Here, Here, according to Eq. (2), the following 
asymptotics arising due to the ion effective mass growth at low frequencies 
holds: 

Im V 1/2 .ON 

cx u; / , (3 

Re V 

As to the measurements reported in Ref. [3] they on the one hand reveal a 
substantial polaronic effect (the ion mass proves to be several times larger 
than rrf'^). On the other hand, at low frequencies Eq. (3) does not hold and 
instead the limiting behaviour 

(Im v/Re v)^^o oc u, (3a) 

is observed which is typical of clusters with fixed mass. An acceptable trade- 
off between (2) and (3) is that the nf^^iuj — >• 0) oc a;~^/^ dependence should 
reach a saturation in the vicinity of ujt ~ 1, i.e. in the frequency range where 
Re F{uj) ~ Im F{ijj). Bearing this in mind, it is easy to show [1] that the 
associated mass m^**(pn, it!„) for the superfluid helium transforms to 

mr(Pn,i?n) = 2.l7ri?^p„, (4) 

which means that the mass becomes frequency independent and its observed 
temperature dependence is due to the factor Pn{T) in Eq. (4) and its maxi- 
mum value substantially exceeds the ideal associated mass m^d- However, the 
assumption Re F{u) ~ Im F{u) is not quite consistent with the conditions 
of the experiments reported in Ref. [3] . 
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Actually the effective mass saturation in the limit rrf^^ioj — > 0) is related 
to the divergency in the linear approximation of the integral 

W = j p{r)u\r)d\ (5) 

where u{r) is the velocity field falling off anomalously slow (as r~^) at large 
distances from the body moving with a constant velocity v through the liquid. 
The details of this non-linear scenario clearly indicated in Ref. [1] have not 
yet been studied. It is the purpose of the present paper to consider this 
scenario and discuss its applicability to data of Ref. [3] . 
1. Divergency of Eq. (5), and hence of the quantity m"^* 

m:r^V2 = W, (5a) 

can be eliminated in the so-called Oseen approximation [1,7,8] which reveals 
that the behaviour u{r) oc following from the linearized Navier-Stokes 
equation and resulting into the divergency (5) is actually replaced by the 
exponential decay of the velocity field at distances r > R/Re, Re = Rv/t] <^ 
1. Because of this exponential decay the quantity m^j'* becomes finite but 
acquires a non-linear dependence on velocity v 

mil' - ^ (l/Re), Re ^Rv/r)<^ 1, (6) 

where rriid is defined in Eq. (2) and Reis the Reynolds number. The additional 
factor in (6) accounts (with the logarithmic accuracy) for the presence of a 
laminar trace generated by a sphere moving through a viscous liquid. 

The result (6) is consistent with both the linear asymptotics rrf^^ioj —>■ 
0) oc a;^^/^ and the requirement of reaching a plateau in the limit. The same 
treatment reveals the physical reasons of the divergency and the mechanism 
of its elimination. However, the progress in understanding leads also to some 
formal "losses" here since the cluster dynamics in the most interesting domain 
mil' /irtifi ^ 1 becomes non-linear. In particular, the Fourier representation 
providing a transparent interpretation of the difference between (3) and (3a) 
can no longer be used. Analysis of dynamic properties of the cluster described 
by Eq. (6) requires an alternative approach. In the present paper we discuss 
with this aim the problem of a step-like external force F acting on a particle 
in viscous liquid. An appropriate quantity for the efficient particle mass is 
the combination 
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Figure 1: Crossover from the A(t) curve defined by the asymptotics (12) 
(thick sohd line) to the long-time behaviour described by Eq. (15). Thin 
solid, dotted, and dashed lines correspond to 6/a = 3, 2 and 1.4, respectively. 



which is calculated below. 
For the Drude dynamics 

{mo + mid){v + v/T) = eE{t), (8) 

the problem of a step-hke external force eE apphed to the particle leads to 
the velocity v{t) 

v{t) ^Voo[l - exp{-t/r)], Voo = eEr/{mo + rriid), 

ADrude = l/r, (9) 

where mo is the cluster bare mass, rriid is defined by Eq. (2), and e,E — const 
are the ion charge and driving electric field. 

On the other hand, the conventional Langevine equation 

(mo + mid)v + 67rRr)v{t) = eE{t) (10) 

yields 

v{t) = foo[l - exp(-i/T*)], = eEn/{mo + rriid), 
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Figure 2: Temperature dependence of the viscosity r) and inverse mobility of 
negative ions /^"^ in helium. 
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Figure 3: Experimental (full squares) and theoretical (solid line calculated 
according to Eq.(17)) effective mass of negative ions in helium 



(mo + rriid) 

^Langev — 1/ ''"* 

In both cases (Eqs. (9) and (10)) the combination A does not depend on 
time, and for the Drude case it also does not depend on the ion mass. 

It is natural to develop the Stokes dynamics which we are interested in 
on the basis of the linear approximation where the Fourier component F{uj) 
of the efficient drag force acting on the sphere moving in a viscous liquid is 
given by Eq. (1). For the velocity v{t) this force yields [8,9] 

(A -1m 7P r ^^P (^1^) _ (^2^) 



q sf - si si S2 
1 /■* exp [slit - r)] - exp [sl{t - r)] 



— [ 

\pK Jo 



Jo 

v(t = 0) = 
where Si and S2 are roots of the equation 

+ ps + q = 0, 



dr] (11) 
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? = , , P = 3— , 

mo + rriid mo + rriid 

eE 

7 = , K, — GnRrj. 

mo + rriid 

In agreement with the initial condition the left-hand side of Eq. (11) turns 
into zero at t — (which is readily verified if one employs the formulae 

Sl + S2 = -p, SiS2 = q). 

In the opposite limiting case (where |si2|^ ^ 1) the asymptotics 

holds. The structure of Eq. (12) reveals that the velocity v{f) approaches 
its asymptotic value v[o6) = '~f/q following a square-root law. It is also 
interesting that this asymptotics is formed exponentially with the typical 
time 

Ti,2 ~ (13) 

The appearance of exponentials (13) has a transparent qualitative inter- 
pretation. At the initial stage of the adjustment of v{t) to its steady-state 
value the process (11) resembles the Langevine scenario (10) with the typical 
exponential relaxation to the stationary behaviour and the relaxation time 
inversely proportional to the constant cluster mass. Later, when the associ- 
ated mass starts to compete with the bare mass (either mo or rriid), a specific 
square-root approach to the stationary regime develops which is absent in 
the traditional dynamics. 

In addition to being of substantial interest in itself, the results (11-13) 
proves to be important for correct formulation of our main problem of the 
velocity relaxation for the ion with mass (6). The point is that this defini- 
tion does not apply at the initial stage of the process where both v{t) and 
the associated mass are growing with time, i.e. (dm/dv) > following the 
scenario correctly described by Eq. (11). On the contrary, in the situation 
described by Eq. (6) one has (dm/dv) < 0. 

Therefore, at some intermediate stage of the relaxation process the quan- 
tity (dm/dv) should have an extremum. At present we are unable to deter- 
mine its position on the time axis in a self-consistent way, for example by 
solving the harmonic problem in the Oseen approximation. Therefore, the 
suggested approximate solution consists of two parts. The initial stage is 



7 



described by linearized dynamics (11-13). Its final part serves as the initial 
condition for dynamics with mass (6). The matching time and the corre- 
sponding ion velocity f* are taken to be ~ ti and ~ v{ti) where t\ is 
the longer one of the two times ri,2 (13). 
The equation of motion to be solved is 

+ RelnT/Re) ^ + ^''^'^''^'^ ^ 
m^^mo + mu Re = Rv/rj (14) 

where v^, and are the matching velocity and time. To find and t^, one 
should solve Eq. (14) for arbitrary values of and and then to study 
the general conditions of the intersection of curves A(t) resulting from Eqs. 
(11) and (14). Then the intersection domain should be used to determine 
the values of v^^ and t^. In our approximate approach we adopt ~ Ti, 
~ v{ti). 

In the most interesting limit Re ^ 1 Eq. (14) can be simplified and 
explicitly integrated to yield 

I 77\ ~ I =t-U, 

— av[t) b — av* 

GttR^ , eER 

a = , 6= . (15) 

rriid rnidV 

It is obvious that under the conditions b — av{oo) — > the process (15) 
approaches the stationary regime with the velocity 

v{oo) = Vmax = eE/QirRr) (16) 

The general behaviour of A(t) for v{t) specified by Eqs. (11) (in the range 
where asymptotics (12) is valid) and (15) is presented in Fig. 1 where the 
choice of parameters and is also illustrated. 

2. The above analysis allows two qualitative conclusions to be made. 
First, in the domain cut^ <S 1 the cluster associated mass indeed has the 
structure (6) with the efficient velocity v = Vmax reached by the ion in a 
single cycle. In these estimates employing the results obtained for a semi- 
infinite time interval the role of uj~^ is played by appropriate finite time 
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interval containing the initial point. Second, in the same frequency range 
Vmax can be estimated using either the data on stationary ion mobility in 
the Stokes form (16), or direct experimental measurements of that mobility 
covering also the transitional (Knudsen-Stokes) domain. 

Let us discuss the data of Ref. [3] within the framework provided by Eqs. 
(14-16). First of all one should estimate the extent to which the inequality 
UJT^, ^ 1 is satisfied assuming that u < 10~^° s and r^, > Ti (13). The 
available experimental data on the helium viscosity and normal component 
density yield 

?7 ~ 2 • 10-' g/(cm s), T, = PnRl/^ri < 10"^^ s, 

UJTs ~ 10"^ < 1. 

Further, the Reynolds number 

Re ~ PuVdRsM, Re < 1 

proves to be dependent on the driving electric field strength whose value was 
not given in Ref. [3]. It is only clear that , the field should be sufficiently 
weak so that the Reynolds number is small (since otherwise it is impossible 
to explain the observed mass enhancement). Under these conditions, bearing 
in mind that 

= Mu^ 27ri?3pn/3 

and assuming that in all measurements the electric filed amplitude Ed was 
kept constant, 

vd — I-iEd, Ed = const 

one obtains 

Mr-M4 = ^=^ = con.t.!^ (17) 
Rg 6Vd o/iEd a* 

where /i is the stationary ion mobility [3]. Thus, in contrast to Eq. (4) the 
arising interpretation of the efficient mass temperature dependence is related 
to the simultaneous effects of both r]{T) and /u(T). The corresponding plots 
together with the data of Ref. [3] are presented in Figs. 2 and 3. 

To sum up, one can say that the complicated temperature dependence 
of the Stokes associated cation mass is a manifestation of a rather general 
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phenomenon inherent to motion of various mesoscopic clusters through vis- 
cous hquid. In the extreme case of small Reynolds numbers their efficient 
mass proves to possess a substantial velocity dependence. The outlined effect 
should be taken into account in the ion dynamics in various electrolytes as 
well as in the calculations of equilibrium and dynamic properties of colloid 
systems, etc. 
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